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The thermodynamic potential is calculated for a uniform superfluid gas of fermi atoms from the 
mean field BCS equations including corrections from induced interactions, Hartree-Fock energies and 
quasiparticle selfenergies. The entropy, specific heat and sound modes are calculated as function 
of temperature, density and interaction strength from the BCS to the unitarity limit and around 
the BCS-BEC crossover. The second sound speed is of particular interest as it is a clear signal of a 
superfluid component and it determines the critical temperature. 



Recent experiments probe systems of fermions near 
Feshbach resonances by expansion [1-5] , collective modes 
[6,7] and RF spectroscopy [8]. Interesting new strongly 
interacting or dense phases of fermions and bosons are 
created, e.g., that associated with the crossover from a 
superfluid or normal Fermi gas to a molecular BEC. The 
experiments provide strong evidence for a superfluid state 
at low temperatures. The next generation of experiments 
will measure sound velocities in these systems. Second 
sound is particularly interesting as it is a clear signal of a 
superfluid component and it also determines the critical 
temperature. 

Ho [9] has studied the thermodynamics in the unitarity 
limit at the Feshbach resonance and expressed a number 
of thermodynamic quantities in terms of unknown uni- 
versal parameters. The purpose of this work is to extend 
the analysis to all temperatures, densities and interac- 
tion strengths in the BCS limit and around the unitarity 
limit, and to calculate the thermodynamic quantities in 
terms the binding energy per particle E/N and the pair- 
ing gap A at zero temperature. These quantities are 
known rather accurately as functions of density and in- 
teraction strength from Monte Carlo [10,11] and exper- 
imental results are compatible. These functions will be 
treated as input into a crossover model that is based on 
the mean field BCS equations of Refs. [12-15], however, 
including important corrections in the crossover model 
from induced interactions, self-energies and Hartree-Fock 
energies. We will then calculate the thermodynamic po- 
tential in the superfluid phase from which the entropy, 
specific heat, first and in particular second sound speeds 
can be calculated at low temperatures in the superfluid 
and normal phases from the dilute BCS limit up to and 
around the unitarity limit. 

The mean field BCS equations of Refs. [12-15] have 
become a standard reference for at least a qualitative de- 
scribing the BCS-BEC crossover as function of density, 
interaction strength and temperature. It describes a sys- 
tem of Fermi atoms of mass m with two spin states in spin 
equilibrium n± = n-\ = n/2 interacting through a s-wave 
scattering length a. The two-body interaction range is 
assumed to be short as compared to \a\ leaving only two 
length scales a and the interparticle spacing (or kp 1 ). 
At zero temperature all the physics depend only on one 
variable which is conveniently chosen as x — l/(akp). It 



varies from — oo in the dilute BCS limit through the uni- 
tarity limit at Feshbach resonance at x — to x — > +oo 
in the molecular BEC limit. 

The mean field BCS gap equation 
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is valid for any coupling strength g = —4nh a/m and can 
thus at least qualitatively describe the smooth crossover 
from a BCS state to a molecular BEC. As usual e k = 
h 2 k 2 /2m - /i A and E k = ^/e\ + A 2 . Note that the 
chemical potential //a does not include Hartree-Fock 
energies and the quasiparticle energy does not contain 
any self-energy. The thermal distribution function is 
fk = (cxp(E k /T) + l)- 1 in units where k B = 1. With 
the equation for number density conservation 
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the gap and chemical potential can be calculated as func- 
tion of density, temperature and interaction strength. At 
zero temperature the last term in Eqs. (1) and (2) vanish 
and the gap and chemical potential are easily calculated 
as function of x = l/(ak F ) (see Refs. [13,15] and Fig. 1). 

In the dilute BCS limit A -C /xa — Ef, where the 
Fermi energy is Ep — h 2 k F /2m, and the Fermi momen- 
tum Tikp is given in terms of the density n = k F /3ir 2 . 
The sum in the gap equation is at zero temperature sim- 
ply (mkF/2n 2 ) \h(kEf /Aq), where k — 8/e 2 . The super- 
fluid gap at zero temperature A = A(T — 0) becomes 



Aq = kEf exp 
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Gorkov, however, found that induced interactions lead 
to a higher order correction such that: a -1 — * a -1 — 
2k p ln(4e)/37r [16]. In the above gap equation, where 
the l.h.s. is proportional to x = l/{akp), this correc- 
tion corresponds to adding or shifting x by the amount 
21n(4e)/37r. As result the Gorkov gap has k = (2/3) 7/3 
in Eq. (3). In the following we shall add this constant 
shift correction not only in the dilute BCS limit but gen- 
erally at all x. Therefore the resulting gap and chemical 
potential as shown in Fig. 1 both deviate from standard 
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results by the above shift in x. The resulting gap is exact 
in the dilute BCS limit and generally in good agreement 
with Monte Carlo calculations [10] up to x^l. 

Neither self-energies nor Hartree-Fock energies are in- 
cluded in the mean field BCS equations. Monte Carlo cal- 
culations provide detailed insight in their contributions 
relative to the effect of pairing because pairing can be 
included or excluded in the trial wave functions. The en- 
ergy per particle E/N — (3/5)Ep(l + (3) is expressed in 
terms of the universal function (3(x) = E int /Ek in [17]. In 
the unitarity limit [3{x = 0) = —0.46 without pairing but 
(3(x = 0) = —0.56 with pairing in the Monte Carlo cal- 
culations of Ref. [10] at zero temperature. Therefore, the 
dominant contribution to the binding energy and chemi- 
cal potential does not come from pairing but from other 
correlations in the wave function that contribute to the 
Hartree-Fock energy. We shall therefore distinguish be- 
tween ha, which includes only the pairing effects, and 
the correct chemical potential 
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which includes both pairing and Hartree-Fock energies. 
Both are shown in Fig. 1. On the BCS side and even 
in the unitarity limit yUA differs from Ef by a small 
amount only. This is compatible with the small de- 
viation between the chemical potentials calculated by 
Monte Carlo with and without pairing. On the BEC 
side x = \/{akp) > the pairing contribution to binding 
energies become increasingly important, and the chemi- 
cal potential rapidly drops toward the molecular binding 
energy /lia — A* — ~ h 2 /2ma 2 in the dilute BEC limit. Re- 
cent experiments on expansion energies [1-5] and collec- 
tive modes in traps [6,7] are compatible with the Monte 
Carlo calculations. 

In both the hydrodynamic limit and for a superfluid 
gas the first sound is given by the adiabatic sound speed 
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where vf = fikp/m, [3' = d/3/dx, etc. In the 
dilut e BCS limit an d at low temperature u\ = 
VFyJ^ + (2/7r)a/cf)/3. In the unitarity limit u\ — 
y/2yi{x = 0)/3m = v F + /?(0))/3 ~ 0.37v F - In 

the dilute BEC limit ui — \J (ir/2)h 2 na m /m 2 , where 
a m ~ 0.6a is the molecular scattering length [18,11]. The 
first sound speed from Monte Carlo calculations [11] is 
plotted in Fig. 1. 

In a normal Fermi liquid the first sound speed is ex- 
pressed in terms of the Landau parameters as: u\ = 
(u|/3)(l + F )/(l + Fi/3), at low temperatures. The 
Fermi liquid theory can be generalized to finite tempera- 
tures T^,0.5Tp and arbitrary relaxation times [19]. The 
effective mass m* = (1 + Fi/3)m may be determined 
from, e.g., the specific heat as shown below, and F can 



then be determined from u\ or equivalently f3(x). In the 
dilute BCS limit m*/m = 1 + [8(7 In 2 - 1) / 15n 2 }a 2 k 2 F . 
In the mean field BCS equation m* = m because selfen- 
ergies are not included. We shall in the following assume 
that selfenergies can be included in the quasiparticle en- 
ergies such that they lead to an effective mass at the 
Fermi surface. Eventually m* will have to be measured 
and/or calculated by Monte Carlo as in the cases of A 
and \x. The Landau parameter Fq varies from zero in the 
BCS limit to —1 in the BEC limit where, however, the 
liquid will only be in the normal state for temperatures 
above the critical temperature. For a BEC this condition 
is T > T c = (n/2C(3/2)) 2 / 3 7r/m ~ 0.218-E F - Fermi liq- 
uid theory also describes collisions and the transition be- 
tween the hydrodynamic (first) sound and the collision- 
less (zero) sound as the collision rate decrease. The Lan- 
dau damping can also be calculated and for — 1 < F < 
the zero sound mode becomes purely imaginary in the 
collisionless limit. 
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FIG. 1. The chemical potential /ia and pairing gap Ao in 
the crossover model, and jj, from Monte Carlo calculations [11] 
(all in units of Ef and at T = 0). The first and second sound 
speed at maximum (T ~ 0.7T C assuming m* = m) are also 
shown in units of vf- 

In the following we shall refer to the mean field BCS 
equations with the above described corrections from in- 
duced interactions, Hartree-Fock energies and effective 
mass as simply the "crossover model". It contains the 
correct pairing and mean field in the dilute BCS limit to 
leading orders, and is also correct in the unitarity limit 
to a good approximation. It docs not solve the full many- 
body problem of strongly interacting Fermi gases but 
should regarded as an approximate model that includes 
some of the most important physics. As will shown in 
the following it has the virtue that it produces definite 
predictions for a number of observables which are about 
to be measured. 

The gap at finite temperature is found by solving the 
gap equation. It is a major simplification that we for 
most purposes treat the pairing gap as being small as 
compared to the Fermi energy for x^,0. The next or- 
der correction to, e.g. the gap itself turns out to be 
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(A/2Ep) 2 , which even in the unitarity limit and at zero 
temperature is less than 10% in both the crossover model 
and Monte Carlo calculations. To the same order in 
A/Ep we can also approximate /ia — Ep in the gap 
equation. 

Inserting the zero temperature gap of Eq. (3) back 
into the gap equation (1), we arrive at the expression for 
the finite temperature gap 
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Because the gap is small, only the quasi particle energies 
around the Fermi surface only contribute to the integrals. 
Therefore, the effective mass in the quasiparticle selfen- 
ergy near the Fermi surface appears in the level density. 
The effective masses should also be included selfconsis- 
tently in Ao. Therefore Ao may deviate from the Gorkov 
gap, and Ao should be taken from experiment or Monte 
Carlo calculations as will be assumed in the following. 

The temperature dependence of the gap resulting from 
Eq. (6) is plotted in Fig. 2. It vanishes above the 
critical temperature, T c = (j/tt)A , and has the form 

for T c - T < T c . Near 
V27rA Texp(-A /T). 



A = 7rT c ^8/7C(3) N /T- T/T c 
zero temperature A = Aq 



The thermodynamic functions can be calculated from 
the thermodynamic potential per volume fl = —P. We 
make the standard assumption that the Hartree-Fock 
terms in the superfluid Q s and normal state fi n ther- 
modynamic potentials are the same. The difference is 
then given in terms of the pairing coupling as 
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The first order temperature correction to the ther- 
modynamic potential in the normal phase is Q n = 
-N(0)ir 2 T 2 /3, where N(0) = m*k F /2TT 2 is the level den- 
sity. It should be noted that the above thermodynamic 
potential is assumed to be valid also in unitarity regime, 
where the Hartree-Fock terms are of order (3Ep. How- 
ever, as discussed above the pairing has relatively small 
effect on the chemical potentials and we may therefore 
expect that the difference between Hartree-Fock terms 
in the superfluid and normal phases also remains small. 

Reinserting the gap of Eq. (6) in the gap equation we 
find for the coupling 
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Note that k now included induced interactions and pos- 
sible effective mass corrections such that the correct A 
of Eq. (3) is reproduced. Inserting the coupling of Eq. 
(8) into the thermodynamic potential (7) and again ex- 
ploiting that A <C ^a, it reduces to 



fl s = -N(0) 
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The crossover model thus arrives at the standard ex- 
pression for Q s and therefore also the standard en- 
tropy S s = —(dQ s /dT)v,ij,, and specific heat C s = 
T{dS s /dT)v,fi in superfluid phase. In the normal 
phase S n = C n = N(0)2n 2 T/3. At low temper- 
atures S s = ^(0) v /27rA3/Tcxp(-A /r) and C s = 
iV(0) v /27rA§/r3exp(-Ao/T). Near T c : S s /S n (T) = 
l-(l+0(l-r/T c ) and C s /C n (T c ) = £-3.77(1 -T/T c ). 
Here f = 1 + 12/7C(3) ~ 2.43 is the superfluid specific 
heat relative to the normal one at T c . 

Finally, we need the superfluid (London) density 
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where n = n s + n n is the total density. At low tem- 
peratures n s /n = 1 — v /27rAo/Tcxp(— A /T) whereas 
n s /n = 2(1 - T/T c ) near T c . 




FIG. 2. Thermodynamic quantities for a superfluid as cal- 
culated in the crossover model vs. T/T c . Shown are the gap 
A(T)/A , entropy S 3 (T)/S n (T), specific heat C S (T) /C„(T)£, 
and superfluid density n s (T)/n. The second sound is plotted 
in units of VF\Jm* JrnT c / 'Ep. 



We now have all the thermodynamic quantities avail- 
able for calculating the second sound speed 
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as shown in Fig. 2. At low temperatures the second 
sound is linear in temperature 
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Around T ~ 0.7T C the second sound speed has a broad 
maximum of 
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Near the critical temperature T c — T <C T c second sound 
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has the characteristic w 2 oc ^/l — T/T c behavior and the 
vanishing point determines T c . 

Ho [9] predicted the scaling behavior of the above ther- 
modynamic quantities in the unitarity limit (x = 0) 
near T c from univ ersal scaling laws. In his analysis the 
u 2 oc yjl — T/T c dependence near T c and x = was 
found but quantitatively it depended on unknown param- 
eters. In comparison the crossover model can be applied 
to all temperatures, densities and interaction strengths as 
long as x^O. The results are given in terms of the func- 
tions Ao(a:), n{x) and m*(x), of which the first two are 
known from Monte Carlo calculations and experiments 
as discussed above. 

The crossover model breaks down in the BEC limit, 
x^0.5, because the gap becomes large, the chemical po- 
tential becomes negative, gaussian fluctuations become 
important at T c [15], and pseudogaps appear [20,21]. The 
second sound with maximum u 2 — 0.53vf(T c /Ef) can 
therefore not be trusted on the BEC side for x^0.5 in 
the present model. If the first and second sound con- 
tinue to decrease and increase respectively as we go from 
the BCS to the BEC limit (see Fig. 1), i.e. as x>0, they 
would cross. However, the two sound modes are coupled 
and undergo avoided crossing. This is known to occur in 
a BEC at the temperature T* — \g\n [22,23]. Below this 
temperature (T;$T*) the Bogoliubov condensate mode 
u 2 = \g\n s (T)/m is first sound whereas above it is sec- 
ond sound. Above t he fir st sound speed is the classical 
result Mi = 0.925-v/T/m. At the crossings, where the 
normal and superfluid components couple strongly, we 
may expect stronger damping of both sound modes. 

In summary, the mean field BCS equations were cor- 
rected with induced interactions, Hartree-Fock energies 
and self-energies. The resulting crossover model is ap- 
proximately correct at all temperatures, densities and 
interaction strengths up to x^,0.5, and exact in the di- 
lute BCS limit. A number of thermodynamic quantities 
were calculated for gas of Fermi atoms with superfluid 
and normal components. Second sound is particularly 
important because it only appears in a superfluid gas 
and it reveals T c where it vanishes. The predicted sound 
speeds in the superfluid phase can be tested in upcoming 
experiments. The validity of crossover model will be put 
to a test in the unitarity limit and measurements of the 
sound modes on the molecular BEC side will provide im- 
portant information on new phenomena as pseudogaps, 
quasiparticle energies, etc. 
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